Superposition states of ultracold bosons in rotating rings with a realistic potential barrier 
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In a recent paper [Phys. Rev. A 82, 063623 (2010)] Hallwood et al. argued that it is feasible to create large 
superposition states with strongly interacting bosons in rotating rings. Here we investigate in detail how the 
superposition states in rotating-ring lattices depend on interaction strength and barrier height. With respect to 
the latter we find a trade-off between large energy gaps and high cat quality. Most importantly, we go beyond 
the (S-function approximation for the barrier potential and show that the energy gap decreases exponentially 
with the number of particles for weak barrier potentials of finite width. These are crucial issues in the design of 
experiments to realize superposition states. 



PACS numbers: 67.85.Hj; 03.75.Gg; 03.75.Lm 

Ultracold bosons in rotating ring traps are a fascinating sub- 
ject of research because the particles in this system can form 
various multi-particle superposition states. Some of these su- 
perposition states may play a central role in studies of the 
quantum-to-classical transition or in applications such as 
entanglement-enhanced metrology ■ 

Most of the work in this area has focused on the produc- 
tion of so-called NOON states in uniform ring lattices 01 
and superlattices ff}. In these systems there are critical rota- 
tion frequencies at which the single-particle spectrum is de- 
generate so that even weak interactions lead to strong coiTela- 
tions between the particles. The energy gap between ground 
and excited states for these systems does, as one might expect 
for multi-particle transitions, decrease exponentially with the 
number of particles. This does, of course, limit the schemes 
proposed to generate these quasi-momentum NOON states to 
operate with a modest number of atoms . 

Very recently, the authors of Ref . showed that for strong 
interactions the ground state of the Lieb-Liniger model with 
a (5-function potential barrier is a superposition of states with 
different total (angular) momenta of the particles. They found 
that in this regime the energy gap is independent of the num- 
ber of particles for small barrier heights and proportional to 
the number of particles for strong baiTier heights. They come 
to the conclusion that the production of superposition states 
with much larger number of particles will become feasible. 

In this paper we study in detail the different superposition 
states that occur for ultracold bosons in rotating ring lattices 
with a realistic potential barrier. We show that these superpo- 
sition states strongly differ in the regimes of weak and strong 
interaction strength (see also Ref. |9]) as well as of weak and 
large potential barriers. Most importantly, we demonstrate 
that for weak barriers the energy gap decreases exponentially 
with the number of particles as soon as the finite width of the 
potential is taken into account. These are important limita- 
tions to producing superposition states that must be taken into 
account in planning experiments. 



We focus on systems at a critical rotation frequency where 
the single-particle spectrum of the uniform ring is doubly de- 
generate and introduce a potential barrier that breaks trans- 
lational invariance leading to superposition states of different 
quasi-momenta. We study systems with much less than one 
atom per site in order to avoid the superfluid-to-Mott-insulator 
transition for strong interactions. At low filling factors the par- 
ticles only feel the quadratic part of the lattice dispersion so 
that our lattice simulations should be relevant for a continuous 
loop as studied in Ref. |9]. 

We use both exact diagonalization of systems with small 
numbers of atoms and sites and the Bose-Fermi mapping in 
the Tonks-Giradeau regime to characterize the many-body 
ground state. Going beyond Ref. f^, we calculate the ground- 
state overlap with three different trail wavefunctions, the spec- 
trum of the single-particle density matrix, and the experimen- 
tally accessible correlation functions, such as momentum dis- 
tribution and momentum noise correlations. 

Like Ref. we find that the superposition states strongly 
depend on the interaction strength. For weak potential barri- 
ers and weak interactions it is a Bose-Einstein condensate, for 
intermediate strength it is a strongly correlated NOON state, 
and for strong interactions, where the particles are locked to- 
gether, it is the superposition of two center-of-mass motional 
states. In the lattice case we consider that the superpositions 
also depend on the potential height. We find that they are 
degraded significantly if the potential barrier is large, which 
leads to a trade-off between the energy gap and cat quality. 

The limit of strong interactions and weak potential barrier is 
most attractive for generating superpositions with many parti- 
cles. In this case, the energy gap is independent of the number 
of particles for a single-site baiTier [9J. However, as soon as 
we consider a physical potential with finite width, the energy 
gap decreases exponentially with increasing number of parti- 
cles. Since any real potential will have a finite width this will 
give a physical limit to the number of particles in the superpo- 
sition state. 
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Hamiltonian. We consider a system of N ultracold bosons 
with mass M confined in a ID ring lattice of L sites with 
lattice constant d. The ring is rotated in its plane with angular 
velocity f2. In the rotating frame the Hamiltonian is l6l ll0ll 
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where fij — a}-dj and dj are the number and bosonic annihi- 
lation operators of a particle at site j, 9 — AlflLd^/h is the 
phase twist induced by rotation, J is the hopping energy be- 
tween nearest-neighbor sites, U the on-site interaction energy, 
and Vj describes the potential barrier at site j. 

To understand the effect of rotation on the atoms in the ring 
lattice, we write the many-body Hamiltonian ([T]i in terms of 
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the quasi-momentum operators bq — L^j=i "-j^ 
where 2'Kq/dL is the quasi-momentum and g = 0, . . . , L — 1 
an integer In this basis the Hamiltonian ^ has the form iHI^ 
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where -Eg — — 2 Jcos(27r(7/L — 0) are the single-particle ener- 
gies, Vqqi is the Fourier transform of the single-particle poten- 
tial Vqq, = i y,e2'''(«~«'^-''/^, and the modulus is taken 
because in collision processes the quasi-momentum is con- 
served up to an integer multiple of the reciprocal lattice vector 
27r/(i, i.e. modulo Umklapp processes. In absence of a poten- 
tial barrier, i.e. Vj = 0, the single-particle spectrum is twofold 
degenerate for 6 — iv/L which we will refer to as the critical 
rotation frequency or critical phase twist. 

Ground state at the critical phase twist. In the following 
we will focus on the ground state at the critical phase twist 
for a small single-site potential barrier (Vj — for j ^ 0, 
i.e. Vqq' = Vb ^ J)- Note that a single-site barrier leads to 
constant momentum transfer in the first Brillouin zone. Simi- 
lar results were found in Ref. |9J which we will extend here. 

In Fig. [T]we show our results for the many-body spectrum 
and detailed characterization of the ground-state wave func- 
tion using exact diagonalization of a small system for iV = 3 
atoms on L = 5 lattice sites. We can clearly distinguish three 
regimes: weak, intermediate, and strong interactions. 

In the non-interacting system, U ~ Q, all bosons occupy 
the lowest-energy single-particle state. In the presence of the 
potential barrier, Vb ^ 0, translation symmetry is broken and 
quasi-momentum ceases to be good quantum number For a 
weak barrier, Vb ^ J, the ground state of the single-particle 
Hamiltonian is an equal superposition of the quasi-momentum 
states \q = Q) and \q = 1), and the many-body ground state is 
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where |0) is the many-body vacuum. This is corroborated 
by our exact-diagonalization results in the weakly-interacting 
regime: the momentum distribution is {fiq) = (b^qbq) = N/2 
for q = {0, 1} and (hq) = otherwise. The distribution 
of total momentum P{K), i.e. the probability that the many- 
body system has total quasi-momentum 2-KK/dL is binomial. 
It can be calculated from P{K) = \K, n) {K, n\, where 
I A', n){K, n\ is the projector on the Fock state with K total 
quasi-momentum quantum number and n is a shorthand for 
the remaining quantum numbers |9]. Finally, the spectrum 
of the single-particle density-matrix (SPDM) {d\dj) has one 
eigenvalue of size N which clearly signals condensation. 

As soon as the interactions overcome the energy splitting 
induced by the potential barrier, the near degenerate many- 
body states are strongly mixed leading to correlations between 
the atoms. We discussed this mechanism in detail in the con- 
text of NOON-state production in rotating ring lattices and 
superlattices [11] as well as quantum vortex nucleation Il2ll . 
Neglecting the effect of the weak potential barrier, the non- 
interacting many-body spectrum at the critical phase twist is 
N + 1-fold degenerate. Most of the degeneracy is lifted by the 
interactions in first order of perturbation theory, but the states 
(6q)^|0) and (&|)^|0) remain degenerate. In the presence 
of the potential barrier these two states are coupled at some 
higher order of perturbation theory (even if the system is non- 
commensurate), the degeneracy is lifted, and the ground state 
is the quasi-momentum NOON-state dgt] 
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While the momentum distribution is almost unchanged with 
respect to the non-interacting case, the momentum noise 
correlations {fiqUq') — {nq){nq>) pick up the large fluctua- 
tions in the NOON-state (|4|i. The distribution of total quasi- 
momentum P{K) is bimodal, and the SPDM has two large 
eigenvalues of size N/2 indicating fragmentation of the con- 
densate. Finally, the overlap with the many-body quasi- 
momentum basis states unambiguously proofs that the ground 
state is indeed the quasi-momentum NOON-state (|4]i. 

In passing we note that the presence of NOON-states can 
be observed in time-of-flight expansion after inducing many- 
body oscillations with a quench in the rotation frequency ■ 
Their frequency is given by twice the energy gap which in the 
regime of intermediate interactions is exponentially small in 
the number of particles |l5t|6!]. 

Let us now discuss the strongly-interacting regime. In ab- 
sence of the potential barrier, i.e. Vb = 0, quasi-momentum 
is a good quantum number. Strong interactions (U/ J — > oo) 
lead to fermionization of the bosons within each subspace of 
total quasi-momentum K. In this regime the many-body en- 
ergy spectrum is equal to the single-particle spectrum (with 
(anti-)periodic boundary conditions for (even) odd N [13]) 
which is degenerate at the critical phase twist. The potential 
barrier then breaks translation symmetry and couples the two 
degenerate ground states and li/j^^^^), which are the 

ground states of the quasi-momentum subspaces with K — 
and K — N . We thus call the ground-state wave function a 
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FIG. 1: (Color online) Energy spectrum and characterization of the ground-state wave function for N = 3 particles on L = 5 sites at the critical 
phase twist with Vo/J = 0.01 as a function of the interaction parameter U/J: (a) many-body energy spectrum, (b) momentum distribution 
{fiq), (c) momentum noise correlations (nqfigt) — {nq){nqi), (d) spectrum of single-particle density matrix {a\aj), (e) distribution of total 
quasi-momentum P{K) (as defined in the text), and (f) fidelity with respect to the states \tpo}, and \ipoo) in Eqs. (jSj, (|4j and l|5j. 



center-of-mass superposition 
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Compared to the regime of intermediate interactions the two 
largest eigenvalues of the SPDM are smaller in the strongly- 
interacting limit pointing to the fact that the center-of-mass 
superposition (|5]l is only a partially-fragmented state. Since 
many quasi-momentum basis states contribute to (|5]l the over- 
lap with a quasi-momentum NOON-state (|4]i is also smaller 

Please note that the distribution of total quasi-momentum 
P{K) is not sensitive to the difference between the states (|4]i 
and (|5]l. This is why going beyond Ref. with a more de- 
tailed analysis in terms of SPDM and overlap is important. 

Effect of large barrier heights. So far we have focused on 
the case of a small single-site potential barrier Vq <C J- Let 
us now study how larger barriers affect the system. 

In Fig.|2]we plot the energy gap AE and the cat quality Q = 
4P(0)F(iV) |9] as a function of the bamer height Vq. We find 
that the gap is small for Vb ^ J, increases sharply around 
Vo = J, and levels off for Vq >• J. As Ref. |9] we see that 
the energy gap is independent of the number of particles in the 
small barrier limit and increases with the number of particles 
in the large barrier limit. However, particularly interesting is 
the dependence of the cat quality Q and the distribution of 
total quasi-momentum P{K) as the barrier height increases. 
Fig. 12] and its insets clearly show that the properties of the 
superposition states are degraded if the barrier height is raised 
to increase the energy gap /S.E. This effect becomes more 
severe as the number of particles increases. 
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FIG. 2: (Color online) Energy gap AE (solid) and cat quality Q — 
4:P{0)P{N) (dashed) versus single-site barrier height Vo for L = 
9 as well as N = 2 (blue), TV = 3 (green) and iV = 4 (red) at 
the critical phase twist in the strongly-interacting limit U/J oo. 
Insets show distribution of total quasi-momentum P{K) for Vo/ J = 
lO-"* and Vo/J = 10+* for TV = 4. 



While we study a lattice system which is different from the 
continuum model considered in Ref. we expect that our 
results qualitatively carry over to the continuum problem as 
they are obtained at densities well below unit filling. 

Number scaling for finite-width potential barrier In fu- 
ture experiments the barrier will most likely be provided by 
a blue-detuned laser beam with Gaussian beam waist. That 
is why we now generalize the case of a single-site poten- 
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FIG. 3: (Color online) Energy gap AE versus particle number A'^ for 
L = 50, Vo/J = 0.01, b = ^/d = (blue), b^^/d^ 2/3 (green), 
b = £,/d — 4/3 (red), and 6 = ^/d — 2 (cyan). Points are the exact 
single-particle spectrum and lines from perturbative expression 
The inset shows the barrier potential Vj for L = 50, Vo/J — 0.01, 
and6 = C/d = 2. 

tial barrier and consider a Gaussian potential barrier of finite 
width ^, i.e. Vj = Af^^Voe^'^ J"/2C" with the normalization 

A/" = X]^=i s ■ The matrix elements Vqq' are given 

by the discrete Fourier transform of the barrier potential Vj 
which for c? ^ ^ becomes the continuous Fourier transform 

Vqq> = yoe-2'^'(«/'^^)'(«~«')V-^ while for ^ ^ we recover 
the single-site barrier limit Vqqi — Vq/L. 

We can get a perturbative expression for the energy gap us- 
ing the Bose-Fermi mapping [13] and degenerate perturba- 
tion theory first order in Vq. The energy gap AE between the 
ground and first excited state in the Tonks limit for N particles 
is 

A^ = 2^e-2"'(«/'^^)'^'. (6) 
L 

We note that in this weak baiTier limit the exponential scaling 
is independent of the barrier height. 

In Fig. [3] we plot the energy gaps AE which open due to 



the potential barrier as a function of number of particles N for 
various baiTier widths ^/d. In the small baiTier Va/J = 0.01 
the accuracy of the perturbative expression (|6]l is excellent. In 
contrast to the single-site potential barrier, i.e. ^/d ^ 1/L, 
a potential baiTier with finite width ^ leads to exponentially 
small energy gaps for increasing number of particles N. The 
physical reason is that we need increasingly large momentum 
transfer to couple the center-of-mass momentum states. Since 
the matrix elements of a finite-width potential barrier decrease 
exponentially for increasing momentum transfer, we recover 
the exponential scaling of the energy gap AE similar of the 
weakly-interacting regime. 

The NIST group ifTill recently carried out experiments in a 
toroidal trap with radius of 20/im and a 4.3/xm (1/e^ radius) 
repulsive barrier This corresponds to £,/dL w 1/60 and a 
suppression of the gap (|6]l by e^*^^/^"' This will severely 
limit cat state production beyond few tens of particles in the 
weak barrier regime where the cat quality is good. 

Summary. In conclusion, we have shown that a weak po- 
tential barrier in a rotating ring lattice loaded with ultracold 
bosons at low filling fraction supports several qualitatively 
very different superposition states of different total quasi- 
momentum: a condensate for weak, a NOON-state for inter- 
mediate, and a center-of-mass superposition for strong inter- 
actions. Raising the potential barrier will increase the energy 
gap but degrade the superposition states. The energy gap in 
the strongly-interacting limit is only independent of the num- 
ber of particles as long as one can neglect the width of the po- 
tential barrier Otherwise one recovers the exponential scaling 
with increasing number of particles. These are crucial issues 
in the production of superposition states in rotating atomic 
rings. They severely limit the size of superposition states in 
any realizable experimental system. 
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